Double-Pionic Fusion of Nuclear Systems and the ABC Effect — 
Approaching a Puzzle by Kinematically Complete Measurements 
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The ABC effect — a puzzling low-mass enhancement in the irn invariant mass spectrum — is 
well-known from inclusive measurements of two-pion production in nuclear fusion reactions. Here 
we report on first exclusive and kinematically complete measurements of the most basic double 
pionic fusion reaction pn — > d7r°7r° at beam energies of 1.03 and 1.35 GeV. The measurements, 
which have been carried out at CELSIUS- WASA, reveal the ABC effect to be a (7T7t)j=l=o channel 
phenomenon associated with both a resonance-like energy dependence in the integral cross section 
and the formation of a AA system in the intermediate state. A corresponding simple s-channel 
resonance ansatz provides a surprisingly good description of the data. 

PACS numbers: 13.75.Cs, 14.20.Gk, 14.40.Aq, 14.20.Pt 



The ABC effect — first observed by Abashian, Booth 
and Crowe [l[ in the double pionic fusion of deuterons and 
protons to 3 Be — stands for an unexpected enhancement 
at low masses in the invariant 7T7t mass spectrum M^. 
Follow-up experiments @, 0, i, B S 0, B S M, ED] re- 
vealed this effect to be of isoscalar nature and to show 
up in cases, when the two-pion production process leads 
to a bound nuclear system. With the exception of low- 
statistics bubble-chamber measurements [4, 8] all exper- 
iments conducted on this issue have been inclusive mea- 
surements carried out preferentially with single- arm mag- 
netic spectrographs for the detection of the fused nuclei. 

Initially the low-mass enhancement had been inter- 
preted as an unusually large 7T7t scattering length and 
evidence for the a meson, respectively [l|. Since the ef- 
fect showed up particularily clear at beam energies cor- 
responding to the excitation of two As in the nuclear 
system, the ABC effect was interpreted later on by a t- 
channcl AA excitation in the course of the reaction pro- 
cess leading to both a low-mass and a high-mass enhance- 
ment in isoscalar M„, spectra [H, Q EE HE E3] • In 



fact, the missing momentum spectra from inclusive mea- 
surements have been in support of such predictions. 

In order to shed more light on this issue and com- 
plementing our previous result on the double-pionic fu- 
sion to 3 Be [TH, H3] exclusive and kinematically overcon- 
strained measurements have been carried out on the most 
basic system for double-pionic fusion, the pn — > d7r°7r° 
reaction. To this end we have measured this reaction in 
the quasi-free mode pd — ► p spectator dir a ir° at beam en- 
ergies T p = 1.03 and 1.35 GeV at CELSIUS using the 
■In WASA detector setup including the pellet target sys- 
tem [l8|]. The latter provides frozen deuterium pellets of 
size ss 20/im, which cross the beam perpendicularly with 
a frequency of ~ 7kHz. The beam energies have been 
chosen to be in the region of the ABC effect as known 
from inclusive measurements. The experimental results 
on the pd ^ 3 Be7r°7r° and pd — > 3 Be7r + 7r~ reactions and 
first results from the measurements of the dd -^ 4 Be7T7r 
reaction are given in Refs. [19l. |20L |2l|. 

The deuterons emerging from the reaction of inter- 
est here have been detected in the forward detector and 
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identified by the AE-E technique using corresponding in- 
formations from quirl and range hodoscope, respectively 
(see, e.g. Fig.l of ref. [HI). Gammas from the tt° decay 
have been detected in the central detector. 

This way the full four-momenta have been measured 
for all particles of an event — except for the very low- 
energetic spectator proton, which did not reach any ac- 
tive detector element. Thus kinematic fits with 3 over- 
constraints could be performed for each event. 

Due to Fermi motion of the nucleons in the target 
deuteron the quasifree reaction proceeds via a range of 
effective collision energies. Based on the measured en- 
ergies in the exit channel and the thus event-by-event 
reconstructed total energies in the pn system the data 
have been binned into small ranges of effective beam en- 
ergy, in order to reduce the kinematical smearing and 
also to allow the extraction of the energy dependence in 
the total cross section. 

The absolute normalisation of these data has been ob- 
tained by relative normalisation to the quasifree single- 
pion production pn — > dir° [22| measured simultaneously 
with the same particle trigger. For the higher beam en- 
ergy we used in addition the pn — > dn production process 
for an alternative calibration. Both calibration methods 
agree within 10 %. However, since at the higher energy 
the phase space coverage is less complete than at the 
lower energy, the determination of the absolute cross sec- 
tion depends somewhat on the model used in the MC sim- 
ulation for acceptance and efficiency correction. Hence 
the cross sections derived for the higher beam energy are 
plotted in Fig. 4 with an increased uncertainty. Also in 
order to avoid contaminations due to the increased 7T7T7r 
production and other background we introduced in the 
analysis of the high energy data kinematic constraints on 
the spectator proton before kinematic fits. 

Results of our measurements are shown in Figs. 1 - 
4. Since the two ir° particles emerging from the reaction 
of interest are identical particles, observables depending 
only on a single ir° like M^o are calculated by averag- 
ing over both possible combinations. Figure 1 displays 
the Dalitz plots of the invariant mass-squares Mj D ver- 
sus M\ and of the invariant mass-square M? of the 
deuteron with one of the pions versus the same quantity 
with the other pion — both for data and Monte Carlo 
(MC) simulations of a model ansatz discussed below. The 
Dalitz plots are far from being flat, i.e. far from being 
phase-space like. They rather exhibit very pronounced 
enhancements in the region of the A resonance and at 
the low-mass kinematic limit of M 2 a . This unusual 
phenomenon is known as the ABC effect. In this reac- 
tion the 7r°7T channel is free of both isovector (J = 1) and 
isotensor (I — 2) contributions due to vanishing isospin 
coupling coefficients. Hence the observation of the ABC 
effect here means that it must be of isoscalar (I = 0) 
character. Figure 2 depicts the Dalitz-plot projections, 
the spectra of invariant mass-squares M% and M? 
for the two bins of effective collision energies T p = 1.00 
- 1.03 GeV and T p = 1.34 - 1.37 GeV. At both energies 
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FIG. 1: Dalitz plot of the invariant mass distributions for 
M^ versus M* 0ir o (top) and M^ versus M^ (bottom) for 
the quasifree reaction process pn — > d7r°7r° in the range 1.00 
- 1.03 GeV (top) and 1.27 - 1-37 GeV (bottom) of effective 
collision energies. Left: data, right: MC simulation of the 
reaction model, see text. The color-coded scale (z-axis) is in 
arbitrary units 



low-mass M^o^o enhancement and A excitation, respec- 
tively, are the predominating structures. 

Angular distributions are shown in Fig. 3. The distri- 
bution of the opening angle between the two pions 
in the rf7r°7r° (or equivalently pn) center-of-mass system 
(cms) peaks at small angles, in particular if we select 
events with M^o n o < 0.32 GeV/c 2 . This means that the 
low-mass enhancement is associated with pions leaving 
the interaction vertex in parallel. The distributions of 
the deuteron polar angles x in the d7r°7r° cms and 

of the pion polar angles 6^°^° and 0^5 in 7r°7r° and dir° 
subsystems, respectively, are anisotropic and essentially 
symmetric about 90°. The anisotropy observed for the 
latter corresponds just to the one expected from A decay. 
The anisotropy in the 7r° angular distribution in the 7T7t 
subsystem signals some d-wave admixture. It vanishes, 
if we consider only data with M n o„o < 0.32 GeV/c 2 , i.e. 
in the region of the low-mass enhancement. From this 
we deduce that the enhancement is of scalar nature, i.e. 
in total of scalar-isoscalar nature. 

The observed low-mass enhancement is consistent with 
the findings in previous inclusive single- arm measure- 
ments, where only the momentum of outgoing fused nu- 
clei was measured. Note that the momentum of the fused 
ejectile is directly related to the invariant 7T7t mass in 
the restricted kinematical range of the measured scat- 
tering angle (see Figs. 3 - 5 in Ref. [l9| and Refs. 

The momentum spectra of inclusive single-arm mea- 
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FIG. 2: Distributions of the invariant masses M^o^o and 
M d7r o from the exclusive measurements (full dots) of the 
pn — > d-K°TT reaction at effective collision energies of 1.00 
- 1.03 (top) and 1.34 - 1.37 GeV (bottom). The shaded areas 
show the pure phase space distributions. Solid and dotted 
curves give AA calculations with and without the assump- 
tion of a s-channel resonance. All curves are normalized to 
the experimental integral cross section. 



FIG. 3: Angular distributions in the reaction pn — > dn n 
at T p =1.03 GeV for the opening angle between the 

two pions, the angle of the deuteron O^" 7T ° - all in the d7r°7r° 

cms - as well as the pion angles 0^o and &% n in d7r° and 
7T7T subsystems (Jackson frame), respectively. The meaning 
of symbols and curves is as in Fig. 2. In addition data (open 
squares) and calculations are plotted also with the constraint 
M nn < 0.32 GeV/c 2 . 



surements [2t [9| exhibit also a strong central maximum 
(corresponding to large M Wff values) , if the deuteron scat- 
tering angle 0^ ab is close to 0°. A possible albeit contro- 
versial explanation has been given in Ref. 1 by associ- 
ating this central maximum with 7/ and iririr production. 
In our Mjttt data, which cover all angles 6^ Qb > 3°, we 
observe no apparent high-mass enhancement. 

The 7T7T low-mass enhancements observed in the exclu- 
sive data for the 7r°7r° channels are much larger than pre- 
dicted in previous A A calculations 0, [H Q [H, [II] . 
As an example we show in Figs. 2-4 (dashed lines) 
calculations in the model ansatz of Ref. [12|], where we 
additionally included the pion angular distribution in A 
decay and the Fermi smearing of the nucleons bound in 
the final nucleus. We have chosen this model ansatz, 
since in addition to its simplicity it gives the smallest 
high-mass enhancement, i.e. is closest to our observation 
in this respect. 

Since on the one hand the available AA calculations 
obviously fail, but on the other hand the data clearly 
show the AA excitation in their Mn^ spectra, a profound 
physics piece appears to be missing in the interpretation. 
As we demonstrate in the following such a missing piece is 
found by inspection of the energy dependence of the total 
cross section displayed in Fig. 4. Shown are the results 
of this work (solid circles) for pn — > d7r°7r° — multiplied 
by the isospin factor of 2 — in comparison with those of 



Refs. [4, 8] for thepn — > dir + ii~ reaction (open symbols). 
The latter reaction is composed of / = and 1 contri- 
butions. The isovector part can be directly derived from 
the known pp — * d7r + 7r° cross sections by use of isospin 
relations 12211. Since the pp — > d7r + 7r° reaction is very well 
described [23| by conventional i-channel calculations for 
the formation of a AA system in the intermediate state, 
also the isoscalar part of this conventional process can 
be derived by applying isospin relations to the interme- 
diate AA system. The thus obtained 1 = part is shown 
in Fig. 4 by the dotted line, whereas the non-interfering 
sum of isoscalar and isovector parts is given by the dashed 
line. Indeed, the 7r + 7r _ data are in reasonable agreement 
with this sum for y/s > 2.5 GeV. At lower energies, how- 
ever, the measured values are much larger than expected 
from the conventional A A process. 

The so far unexplained structure observed for i/s < 2.5 
GeV in our data is much larger and narrower than ex- 
pected from the conventional AA excitation by t-channel 
meson exchange [l2|, [l3|, [ll| . The cross section maximum 
is shifted substantially below the nominal AA threshold 
being, however, still 200 MeV above the d7T7r threshold. 

Since isospin is very unlikely to be broken on such 
a large scale, a mechanism different from the t-channel 
AA excitation must be the reason for this resonance- 
like structure in the isoscalar sector. These consider- 
ations led us to the concept of a s-channel resonance, 
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FIG. 4: Energy dependence of the total cross section for the 
pn — > dn + -K~ reaction from threshold (y^s = 2.15 GeV) up to 
y/s = 3.5 GeV. Experimental data are from Refs. [§] (open 
squares) and Q] (open triangles) . The results of this work for 
the 7r°7r° channel — scaled by the isospin factor of two — are 
given by the full circles. Dashed and dotted lines represent 
the cross sections for tt + tt~ and channels, respectively, 

as expected from the isovector tt + tv° data by isospin relations 
(see text). The solid curve includes a s-channel resonance in 
the AA system adjusted to describe the ABC effect in the 
7r°7r° channel. 

which couples to isoscalar pn and AA configurations. 
We note that such a kind of dibaryon resonance has 
been predicted by a number of theoretical investigations 

H ill MM m HI- 

For the description of the observables by such a reso- 
nance we use the following Brcit-Wigner ansatz for the 
resonance amplitude 
A R ~ F(q AA ) * D R * D Al * D A2 (1) 
with mass and width of the s-channel resonance being 



m R w 2.36 GeV/c 2 and T R w 80 MeV and where D R and 
.Da stand for s-channel resonance and A propagators, 
respectively. The form factor F(qaa) of the A A vertex, 
which is chosen to be of monopole type, depends on the 
relative momentum (jaa between the two As. Since qaa 
— q-rnr, when neglecting the Fermi motion of the nucleons, 
this form factor is reflected directly in the M v „ spectra 
and causes there the ABC effect by suppression of the 
high-mass region. Fitting the cutoff parameter Aaa of 
this monopole formfactor to the data in the M OT spectra 
results in Aaa ~ 0.2 GeV/c which corresponds to the 
size of a deuteron-like object. 

While other theoretical explanations cannot be ruled 
out at this time, with such a simple ansatz we obtain 
a surprisingly good description of the data (MC Dalitz 
plot in Fig.l and solid lines in Figs. 2-4) both in their 
energy dependence and in their differential behavior. In 
previous attempts (lfjl . [301 ] of data interpretations we had 
to impose a strong attraction or even a boundstate condi- 
tion between the two As, in order to obtain a satisfactory 
description of data on the ABC effect. Having learned 
now from the behavior of the total cross section, that 
possibly a s-channel resonance could be the cause, these 
impositions would find a natural explanation. 

We finally note that after the move of the WASA de- 
tector to COSY new data have been taken there very re- 
cently on this issue with two orders of magnitude higher 
statistics. It is expected that these will be able to finally 
solve the long-standing ABC puzzle. 
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